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The partial success of the block renormalization group techniques is analysed in terms of a func- 
tional operator which formalizes the idea of self-replicability of a system in terms of smaller blocks 
which are similar to the original. The mathematical properties of the fixed points of this transfor- 
mation are analyzed. 

PACS numbers: 



I. INTRODUCTION 

The renormalization group (RG) is one of the most 
relevant theoretical tools in many branches of physics. 
Its main idea is that studying the changes of behaviour of 
a system under a scale transformation can provide very 
useful information. Kadanoff gave the most successful 
mental image when he developed the idea of block spins^. 
Consider a 2D lattice of interacting spins, and split it into 
2x2 blocks. Each block, under some conditions, behaves 
like a single spin, only varying the coupling constants 
which quantify the interaction. 

In this spirit, the block renormalization group (BRG) 
starts its procedure by finding the ground state of a 
hamiltonian on a small system. Then, various systems of 
the same type are put into contact, making up a block, 
and the ground state of the same hamiltonian on the 
global system is searched variationally using as building 
bricks the ground states for each part. Then, the full 
procedure is iterated until we obtain the ground state 
for a system of the desired size. Changing the point of 
view, we may say that the ground state of a big system is 
approximated, variationally using as building bricks the 
ground states for each of its parts. 

The BRG has met variable success in practice. The 
practicioners abandoned it in favour of other meth- 
ods, such as the density matrix renormalization group 
(DMRG), a RG technique of high accuracy, but which 
lacks some important features of the RG, such as 
the intimate relation between fixed points and critical 
systems^. These conceptual problems encouraged a re- 
turn to the analysis of the BRG, in order to understand 
the reason of its failures and successes. This study led 
to the development of the correlated blocks renormaliza- 
tion group (CBRG) by Martm-Delgado et al., which was 
rather successful for problems in quantum mechanics in 
ID and 2D, but for unclear reasons^. 

In this work we introduce a functional transformation 
which explains the outcome of a BRG prescription in the 
case of single-body quantum mechanics (i.e.: obtaining 
the ground state of a hamiltonian acting on square inte- 
grable functions on a subset of R"), which we call the 
replica transformation. In intervals of R, it acts follow- 



ing these steps. First of all, we generate two scaled-down 
copies of the original function, and place them on the 
left and right halves of the interval. Now we find the 
best approximation to the original function within the 
subspacc spanned by these two. Once normalized, this 
best approximation is the replicated function, and the 
scalar product with the original will be called the self- 
replicability of the function. 

Once this operation is generalized to sets of functions, 
we show that a BRG prescription, such as the CBRG, 
is successful whenever the low energy spectrum of the 
hamiltonian on each block constitute an approximately 
self-replicable set. The CBRG attained success by play- 
ing with the boundary conditions between the blocks in 
such a way that this requirement was fulfilled. 

This article is organized as follows. Section |TT] intro- 
duces our model problem, i.e.: to obtain the low energy 
spectrum of the hamiltonian of a free particle in an in- 
terval, along with a naive BRG analysis and the explana- 
tion of its failure. In section [TTTI the CBRG prescription 
is reviewed. Section IIVI introduces formally the replica 
transformation and the self-replicability parameter, both 
for individual functions and for sets of them. In sectionfVl 
we prove some rigorous results about self-replicable sets 
of functions. We conclude in sectionlVH summarizing the 
results and providing some hints for future developments, 
such as the extension to many-body problems. 



II. MODEL PROBLEM AND BRG APPROACH 

A. Model problem 

Our model problem was originally formulated by K.G. 
Wilson as a toy model in order to dilucidate the reason 
of the failures of the BRG^. We are asked to obtain 
the low-energy spectrum of the free hamiltonian for a 
spinlcss particle in a ID box. In mathematical terms, the 
lowest cigenstates of the laplacian on an interval. The 
configuration space is discretized into a graph. Then, 
the hamiltonian becomes a N x N matrix related to the 
discrete laplacian on it: = Hi-\,i = — 1, His = 2, 

for i > 1 and i < N. The boundary conditions (be) 
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are specially important. The discrete analogue of the 
fixed boundary conditions is equivalent to setting if 1,1 = 
Hn.n = 2, while for free boundary conditions we have 
if 1,1 = Hn,n = 1. 

The problem can be generalized to an arbitrary graph 
Q. Then, the hamiltonian of a particle with free boundary 
conditions is just Hf ree = D — A, where D is the diagonal 
matrix in which the entry of each vertex is its degree, 
and A is the adjacency matrix, i.e.: Ay = 1 if there is 
a link between sites i and j, and zero otherwise. Fixed 
boundary conditions have a less natural generalization. 
If the graph vertices have uniform bulk degree d, then 
if fixed = dl — A. A potential energy can be included 
by adding some Vi to the diagonal elements. In absence 
of such a potential, both Hf ree and ff/ixed are positive 
defined for any graph. It can be easily proved that Hj ree 
has a homogeneous zero mode for every possible graph. 

Of course, there are more physical problems which lead 
to the same mathematical formulation, e.g.: a vibrating 
string or a tightly bound electron in a lattice. 

B. BRG approach to the problem 




FIG. 1: Above: the exact ground states of the TV = 20 sys- 
tems. Below: the best approximation within the BRG Ansatz. 

experiment with free be, the result is completely satis- 
factory, but in a trivial way. Both ip L and tp are ho- 
mogeneous, and so is the global ground state. Therefore, 
the first lesson to be obtained is that boundary conditions 
may be determinative for the failure of a RG-prescription. 



Let us consider two ID lattice segments, of TV sites 
each. The ground state of the fixed be hamiltonian is 
known for each of them. Now we attempt to obtain 
a variational estimate of the ground state of the com- 
pound segment using arbitrary linear combinations of the 
ground states of each block as Ansatz. 

In more explicit terms, let i/jq be the exact ground 
state for TV sites. Now we define ip L and tp R to be the 
natural extensions to the lattice of 27V sites: tpf = (i/jq ) j 
if i < TV, and zero otherwise, while ip^_ N = {'4 f o)i if 
i < TV, and zero otherwise also. Now we build a vari- 
ational Ansatz, ip = a^ip 1 " + ctRip R and get an effec- 
tive hamiltonian for these two states. Let H t be the full 
hamiltonian for the composite lattice. Now, using Dirac's 
bra-ket notation, we get 
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The problem, therefore, reduces to the diagonalization 
of this 2x2 effective hamiltonian matrix, if e //- Vari- 
ational approaches are always highly dependent on the 
quality of the Ansatz. In this case, the results prove 
it to be surprisingly inadequate. The BRG approxi- 
mation to the ground state energy of the 40 sites lat- 
tice is -E20+20 ~ 0.202, while exact diagonalization gives 

E$ « 0.00587. This means an error 400%. The cause 
of the failure is apparent when we plot, in figure [l] the 
wavefunctions of the exact ground state for TV = 40 sites 
and the best approximation within the BRG Ansatz. 

The boundary conditions force the wavefunctions to 
take the value zero at the borders of each block, thus 
making a spurious kink appear in the center of the com- 
plete system. On the other hand, if we make the same 



III. REVIEW OF THE CORRELATED BLOCKS 

RG 

A successful real space prescription for this problem 
which respected the BRG spirit was given by Martfn- 
Delgado et al4£. We will describe briefly the method in 
this section. A thorough explanation of the method can 
be found in the PhD dissertation of one of the present 
authors^. 

Let us consider a linear chain of TV sites. We will split 
this chain into TV/, = 2 k blocks of m sites each, for some 
integer k. Each block is labeled by the index p £ [L.iVj,], 
isolated from the others and given free boundary condi- 
tions. It will have, therefore, a certain self-interaction 
hamiltonian A p . We obtain the low energy eigenstates 
of these self-interaction hamiltonians and, with them, we 
build a chain of Ansatze, which we will describe hence- 
forth. 

First of all, we write down the effective hamiltonian 
for two neighbouring blocks, using the 2m states as vari- 
ational bricks. The global be will again be set to be 
free, but in order to reproduce correctly the link between 
them we need two types of operators. Following Martfn- 
Delgado et ali^, we define influence operators as those 
which restore the correct boundary conditions between 
the two blocks and interaction operators as those which 
take into account the dynamical aspect of the joining. 

This superblock hamiltonian, or level- 2 block, is now 
diagonalized exactly, and we retain, out of the 2m eigen- 
states, the m with lowest energy. Those states will rep- 
resent the block as we proceed to the next RG step. Now 
two such level-2 blocks are put together, and the same 
process is repeated. The RG iteration continues until the 
full system is contained in a single superblock. 
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The same technique works, with very reasonable suc- 
cess, with fixed boundary conditions, in presence of a po- 
tential or, with the necessary generalizations, in the 2D 
case. A reference which describes the technical details is 
the PhD dissertation of one of the authors^. 



IV. SELF-REPLIC ABILITY 

The main thesis of this work is the following: the 
reason for the success of the CBRG prescription, as de- 
scribed in the previous section, is the special self-similar 
properties of the eigcnfunctions of the hamiltonian when 
free boundary conditions are applied. 

Given a (normalized) function <f> e L 2 [0, 1], let us de- 
fine the operators L and R as: 



where the matrices a and (3 are given by the scalar 
products 



(5) 



We should emphasize that the replica transformation 
acts on functional subspaces. The considered initial set of 
m functions is just a basis for the relevant subspace, and 
equation [3] for the best approximation is correct only if 
the basis is orthonormal. The scalar products among the 
elements of the replicated set 1Zq4> are given by 



Lcj)(x) = 



V2(f)(2x) if x < 1/2 
if x > 1/2 



R(j>{x) 




1) 



if x < 1/2 
if a; > 1/2 



i.e.: they return reduced copies which are similar to the 
original function for each part (left and right), and the 
\/2 factor is included so that Lcj) and R<j> have the same 
L 2 norm as <f>. We may try to reproduce the original 
4> within the subspace spanned by Lcj) and Rcj>. Is this 
possible? 

Let us consider all functions to be L -normalized and 
let us define the replica transformation: 



Let us remind that J2k l-^fc) (L(j>k\ is a projector on 
the subspace of the left copies, so we define projector 
Plr on the full subspace spanned by the set {Lfa, R<f>i}, 
and observe that 



(n &\n <p 3 ) = (<j>i\p L R\<l>j) 



(6) 



Therefore, the replicated functions are orthogonal if 
and only if the operator Pl,r acts trivially on the original 
set, i.e.: if the junctions are self-replicating. Otherwise, 
we must apply a Gram-Schmidt procedure, let us call it 
G. This way, the wc define the full replica transformation 
as: 



1Zcj> — Af[ {(j)\L<t>) Lcj) + {cp\Rcj>) Rcj> ; 



(2) 



where Af is a normalization constant. Since (L(f>\R(f>) = 
0, this is the best approximation to the original function 
within the given subspace. Its accuracy shall be given by 
the parameter 



s = (R4\<f>) 



(3) 



where the symbol S stands for self-replicability. The 
value 1 means perfect, and means that the two states 
are orthogonal. 

The procedure is easily extended to sets of functions 
or, better, to functional subspaces. Let us denote any 
such set, with m functions, as 4> = {4>i}YLi- which we 
will assume to constitute an orthonormal set in L 2 [0, 1]. 
These functions are approximated within the subspace 
spanned by the 2m functions {L<pi, R(f>i}'^l 1 , which also 
make up an orthonormal set. All the 2m functions, L<p 
and R(j) contribute to the reconstruction of the parent 
wavefunctions <j). We can give a preliminary definition of 
the replica transformation as 



K (, 
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(otijL(j)j + pijR<pj) 



(4) 



n<t> = gk Q (, 



(7) 



In order to perform numerical experiments, we should 
also define the transformation for discretized functions. 
We will assume that 1Z acts internally on a certain dis- 
crete functional space, isomorphous to R*. Therefore, 
when applying the L and R operators, two values of the 
old function must fit into a single new value. The most 
symmetric way of doing this is the local averaging: 




if i < N/2 
otherwise 



(8) 



along with an equivalent formula for the right side. 

The generalization to higher dimensions does not pose 
any theoretical difficulty. In 2D, e.g., functions are de- 
fined in the unit square, which is divided into four re- 
gions. Each function gives rise to four children, out of 
which wc shall attempt its reconstruction. 

Another possible generalization is the choice of a 
Hilbert space different from L 2 . By choosing a Sobolev 
space, we ensure that the replicated function and the 
original one are not only similar in value, but also in 
derivatives. This can be an important issue, as will be 
seen in the next section. 
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FIG. 2: The ground state of a ID laplacian with fixed b.c. is 
not self-replicable. 




FIG. 3: The procedure is iterated. The last box represents 
the fifteenth iteration. 



A. Numerical Experiments in ID 

Let us apply the replica operator to the ground state 
of the laplacian with fixed be on a ID lattice. The best 
approximation is given in figure [2j Here S takes the 
value 0.8488, which is not too bad (with a Sobolev-type 
norm, it would be much lower), but it gets worse as the 
replica operator is iterated. Figurc[3]shows us the second, 
third and fifteenth iterations. After some iterations the 
finite resolution of the computer yields a quasi-constant 
function as approximation. This function is exactly self- 
replicable. 

As another example, let us consider the low energy 
spectrum of the ID laplacian with free be. It is self- 
replicable to a reasonable approximation, although not 
exactly, as shown in figure |4] 

The S parameters are 1, 0.9996, I and 0.9957. This 
means that the ground state (flat) and the third states 
are exactly reproduced. Analysis of the weights shows 




FIG. 4: A single step of the replica transformation on the 
lowest energy states of the free be laplacian. 




FIG. 5: The same functions as in figure [4] after 15 iterations 
(once the fixed point has been reached). 

that: 

• The first state is absolutely self-replicable by itself. 

• The second consist of two copies of itself, the first one 

raised and the second one lowered, using the first 
state to this purpose. The finite slope at the origin 
is not correctly represented. 

• The third state only requires two copies of the second 

one. 

• The fourth state is even more interesting. Both the left 

and the right parts are a combination of the second 
and third states. The finite slope at the origin is 
again incorrectly represented. 

The procedure may be easily iterated without excessive 
distortion. The results are shown in figure [5l These 
functions have a rough look, but factors S are not too 
different: 1, 0.9996, 0.9996 and 0.9949. 

It is also illuminating to perform the same experiment 
on the four lowest energy states of the fixed be lapla- 
cian (see figures [5] and [7]) . The values of the S parame- 
ters at the first iteration are not excessively bad: 0.9537, 
1, 0.9452 and 1 (to four digits). But the fixed point 
yields very different numbers: 0.8431, 0.8724, 0.8516 and 
0.8020. 

There is a great wealth of fixed points of this trans- 
formation, but most of them are non-smooth and have 
non-trivial fractal properties. Numerical experiments 



FIG. 6: The lowest energy states of the fixed b.c. laplacian 
after one iteration. 




FIG. 7: The same functions as figure [6j once the fixed point 
has been reached (15 iterations). 

led immediately to a smooth family of exactly self- 
replicating functions: the polynomials. If we apply the 
Gram-Schmidt procedure on {1 }, it is eas- 

ily checked numerically that the self-replicability param- 
eter is exactly one for all the functions. The meaning of 
this result will be explained in detail in section fVl 

B. Numerical Experiments in 2D 

As it was previously stated, the process may be eas- 
ily generalized to 2D, if instead of splitting the interval 
into two regions we part it into four. The case of the 
eigenfunctions of the free b.c. laplacian yields a fixed 
point which is much smoother than in the ID case, as 
it is shown in figure [5] In the fixed be case, we obtain 
a rather different fixed point, as it is shown in figure [9] 
Among the fixed points we have found a great richness 
of structures. Figure [TOl shows a familiar pattern. 

V. ANALYTICAL FUNCTIONS AND 
SELF-REPLICABILITY 

A very simple argument shows that the only self- 
replicable function on an interval which is analytical is 
the uniform function. The replicated function can be 
always written as 

Tl(j>{x) = ax [0 ,i/2] {x)4>{2x) + 0X[i/2,i] (x)4>(2x - 1) (9) 



FIG. 8: Above, wave-functions for the 2D laplacian with free 
b.c. for a 32 x 32 system. Below, the fixed point we reach. 
These last functions are also smooth. 




FIG. 9: Same concept as in figure [8] but for initial wave- 
functions with fixed be. The fixed point is thoroughly differ- 
ent. 



where X[a,b]( x ) ' 1S the characteristic function for the in- 
terval [a, £>]. 

The self-replicability condition is <fi(x) =1Z(j>(x). Then, 
obviously, (f)(0) = 71(f)(0), from where we deduce that 
a = 1. Now we force the equality of all the deriva- 
tives at that point: (f>^(x) = TL<j) (k) {2x) = 2 k (f>^(2x) in 
[0, 1/2]. Now, picking up the point x = again, we get 
<^ fc )(0) = 2 fc (/>( fe )(0), which proves that all the derivatives 
of the function are zero at the origin so, if the function 
is analytical, it must be a constant. 




FIG. 10: 16 x 16 Wave-functions of higher energy than the 
previous ones, which yield a repeated pattern based on the 
Pascal or Sierpinski triangle. It must be remarked that this 
curious structure is not stable, and that it is due to a slight 
asymmetry of numerical origin in the initial states. 
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A more interesting result is obtained when we analyze 
the self-replicability of a set of functions (equivalently, of 
a subspacc). According to equation [7], we can extend 
the previous argument in the following way. Restricting 
ourselves to the left interval [0, 1/2], the self-replicability 
condition reads 

'£a ij $ j {2x) = $ i (x) (10) 
j 

Let us assume that the functions fa; are analytical. 
Then, derivating k times that expression with respect to 
x we obtain 

J^ ai ^f\2 X ) = ^f\x) (11) 
j 

Restricting ourselves to the point x = 0, that equation 
reads 

^a y f (0)4$f»(0) (12) 

3 

Let us assume that (0) ^ for all k. Then, equa- 
tion |12| implies that matrix otij, which is finite, has an 
infinite set of eigenvalues, l/2 fc for all positive k. Since 
that can not happen, we have proved that, if k > m, the 
derivatives must vanish. In other terms, all the fa must 
be polynomials. 

It easy to prove that the subspace spanned by 
{1, x, x 2 , ■ ■ ■ , x m } is self-replicating. Both results to- 
gether allow us to state the following theorem: the 
only analytical family of self-replicating functions are the 
polynomials. 

VI. CONCLUSIONS 

Let us return to the original question: why do the 
free be CBRG work, while if we try fixed be the failure 
is complete? The reason is that the eigenfunctions of 



the free be hamiltonian make up an approximately self- 
replicable set, because they resemble the polynomials. 

So, the success of a BRG approach requires the build- 
ing bricks to be suitable for the problem at hand. The 
boundary conditions, being the main freedom of the RG 
practicioner, should be chosen with care. But the main 
guide should be this: make the block eigenfunctions as 
close to self-replicability as possible. 

To summarize, we have introduced the replica trans- 
formation on a functional space, which attempts to re- 
produce a function or set of functions by the best ap- 
proximation attainable with reduced and scaled copies 
of itself. We have defined the self-replicability parame- 
ter as a measure of the failure of a function to replicate 
itself. Some numerical experiments have shown the com- 
plex structure of the set of fixed points of this transfor- 
mation, although a proof has been provided that most of 
them are non-smooth functions: only the polynomials, 
up to any order n, constitute an analytical self-replicable 
set. 

The idea of self-replicability was used to explain the 
success of the CBRG, because it uses free boundary con- 
ditions in order to split the blocks. This choice is appro- 
priate because the eigenfunctions of the laplacian with 
those be are approximately sclf-rcplicating. 

A very interesting line of further work would be to 
study the possible extension of these ideas to other prob- 
lems where real space RG has been applied, such as 
many-body hamiltonians. The splitting of the system 
into blocks should be done in such a way that the result- 
ing eigenfunctions of the block hamiltonian are approx- 
imately self-replicable. This poses an interesting chal- 
lenge with immediate applicability to the development 
of numerical methods in condensed matter and particle 
physics. 
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